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Abstract. Partial wave scattering amplitudes in baryon-baryon and baryon-antibaryon collisions and am- 
plitudes for the production and decay of baryon resonances are constructed in the framework of the 
spin-momentum operator expansion method. The approach is relativistically invariant and it allows us to 
perform combined analyses of different reactions imposing analyticity and unitarity directly. The role of 
final state interactions (triangle and box diagrams) is discussed. 
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1 Introduction 

To understand strong interactions at low and intermedi- 
ate energies is one of the important tasks when quantum 
chromodynamics is being studied. At large momentum 
transfer, QCD can be used efficiently due to the smallness 
of the strong interaction coupling constant [1]; the low- 
energy domain can be treated using effective field theories 
[2]. The resonance region is much more difficult to access. 
Lattice gauge calculations are capable to reproduce the 
masses of ground state hadrons [3] but excited states and 
their decay properties are difficult to extract from lat- 
tice data. For further progress, systematic experimental 
information seems to be mandatory to identify the lead- 
ing mechanisms responsible for the mass spectrum and for 
the decay amplitudes of strongly interacting particles. 

Recently, considerable progress has been achieved in 
meson spectroscopy, even though a commonly agreed pic- 
ture has not yet emerged. Recent reviews emphasizing dif- 
ferent views can be found in [4-7]. The main sources of 
recent progress were the study of reactions with multi- 
particle final states. The analysis of data on proton-anti- 
proton annihilation at rest resulted in the discovery of 
a number of particles in the region 1300-1800 MeV [8]- 
[14]; the investigation of the proton-antiproton annihila- 
tion in flight led to a large set of new states over the region 
1800-2500 MeV [15]-[18]. It appeared that the majority of 
the newly discovered states are lying on linear trajecto- 
ries against radial excitation number [19]. Such a pattern 
was not predicted by the classical quark model of Godfrey 
and Isgur [20] using a linear confinement potential and ad- 
ditional interactions due to effective one-gluon exchange 
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forces. More recent calculations based on instanton-indu- 
ced interactions [21,22] can, however, be tuned (by choos- 
ing an appropriate Dirac structure of the confinement 
potential) to reproduce the observed mass pattern very 
well. The pattern can be understood, too, within a 5- 
dimcnsional theory holographically dual to QCD (AdS/ 
QCD) [23] which predicts masses to proportional to (N + 
L) where L is the intrinsic orbital angular momentum 
between quark and antiquark and N a radial quantum 
number. In the light-quark meson spectrum, practically 
all expected states are observed. However, there are a few 
additional states which do not belong to these trajecto- 
ries. These states are candidates to be of exotic nature, 
e.g., they could be glueballs or hybrids. 

The situation in the baryon sector is in some sense 
reverse: for baryons, the quark model predicts a much 
larger number of states than that observed experimen- 
tally. So far, the pattern seems to suggest that not all 
degrees of freedom in the three-quark system are realised 
in the spectrum of excited states. Instead, the pattern of 
excited states follows the same (L + N) pattern [24] which 
is observed for mesons. If this is the case, the fact would 
be an important phenomenon in the physics of highly ex- 
cited states. Still, a detailed verification of this statement 
is needed. On the other hand, the main information on 
baryon resonances has come from the ttN elastic scat- 
tering, and one may hope that many new states will be 
discovered in (i) reactions involving strangeness in two 
hadron final states and in (ii) inelastic reactions induced 
by photons or protons with three or four particles in the 
final states (for example, two pion photoproduction). 

The search for new baryon resonances is of topical 
interest and several experiments like CB-ELSA, CLAS, 
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GRAAL, SAPHIR, and SPRING-8 pursue active searches 
using photoproduction as a tool [25-52]. A few new res- 
onances were suggested [53-55] in fits to these data sets. 
Proton-proton collision experiments can provide an im- 
portant source of information on baryon resonances in- 
cluding exotic states (e.g. pentaquarks [56-63]). COSY at 
the Research Center Julich is providing a wealth of data 
on meson production in proton-proton inelastic scattering 
[64-115]. The experiments Anke and COSY-11 covered 
mainly the threshold region, while TOF covers the full dy- 
namical range. At present, the upgraded WASA detector 
is installed at COSY and will provide high-statistics data 
on the production of neutral mesons in NN interactions 
[116]. The data provide stringent information on nucleon- 
nucleon-meson vertices and on the formation of baryon 
resonances. Selected papers can be found in [117-143]. 

The partial wave analysis of such processes cannot be 
carried out without taking into account the final state in- 
teraction. In many processes, the inclusion of the proton- 
proton interaction dramatically changes the description of 
the data [144]. However, a number of important problems 
for such analyses has not comprehensively developed yet. 
Among these problems are a correct treatment of rela- 
tivistic effects and of the contributions of triangle or box 
diagrams. 

In this paper, we present a relativistically invariant 
approach for the partial wave analysis of proton-proton 
interactions. The method is based on the spin-momentum 
operator expansion suggested in [145-148] . The contribu- 
tion of triangle and box diagrams to the meson production 
processes is discussed and certain examples are considered. 

In Section 2, we present the partial wave expansion 
for baryon-baryon and baryon-antibaryon scattering am- 
plitudes. In Section 3, the unitarity condition for fermion- 
fermion partial wave amplitudes is discussed. The angular 
momentum and spin operators for nucleon-nucleon scat- 
tering are introduced in Section 4, the nucleon-nucleon 
partial wave amplitude is constructed in Section 5. In this 
Section, fermion-fcrmion one-loop diagrams and the cross 
section for the two-fermion scattering are calculated. The 
operators for N A production are constructed in the Sec- 
tion 6. Some examples of amplitudes with multi-particle 
final states are given in Section 7. Properties of the trian- 
gle and box diagrams are shortly discussed in Sections 8 
and 9. 



momentum J (we define t — q 2 = (p'j — Pi) 2 )'- 

M(s,t,u)= £ (Wi)Q S ^A9MPi)) 

S,L,L f ,J 
Ml ■HJ 

x (^)oE,W(»)) 4 S ' L%J V)- (i) 

Here, Q is the i-channel operator (the four-component spi- 
nors ip(p) are given in Appendix 1) and (ii, jjb 2 ■ ■ ■ Hj-i, Hj 
are the indices of the rank J operator. 

Another representation is related to the s-channel (we 
define s — (pi + p 2 ) 2 ): 

M( s ,t,u)= J2 (Wi)Q s ^.iA*±W(-p'2)) 

S,L,L',J 

x (r(-P2)Q°t J .»A k MPi)) 4 S ' L ' L ' J) (s) • (2) 

Here, ip c (—p) are charge conjugated four-component spi- 
nors (see Appendix 1) and <3^f J Al/ are the s-channel op- 
erators, where S, L, J are, correspondingly, spin, angular 
momentum and total momentum of the partial wave in 
the s-channel. The notations of momenta are as follows: 

P =pi+p 2 =p'i+p' 2 , k = -{pi-pi), 

X _ PyP^i IP 7 T _|_ /q\ 

The representation (1) is suitable to consider the t- 
channel meson or Reggeon exchanges, while Eq. (2) is 
convenient for the s-channel partial-wave analysis. The 
representations (1) and (2) are related to each other by 
the Fierz transformation [149], with a corresponding re- 
expansion of the spin-momentum operators. 

In terms of the SLJ representations, the states are 
usually described as 2S+1 Lj. The P-parity can be cal- 
culated as P = (-l) i+1 and C = (-1) L+S . The states 
with S = arc unambiguously defined and they form a 
set of states with J PC = _+ , 1+", 2"+ . . . The states 
with S = 1 and L = J arc also uniquely defined and 
form the set J PC = 1++, 2 — , 3++ . . . The states with 
5 = 1 and L = J — 1 and L = J + 1 have the same 
J PC and can mix with each other, that arc the states 
J PC = 0++,1 — ,2++,... 



2 Selection rules for baryon-antibaryon and 
baryon-baryon scattering amplitudes 

2.1 Baryon pairs with isospin I — 

First, consider the baryon-antibaryon scattering ampli- 
tude in a isospin singlet configuration, for example, the AA 
scattering amplitude. One can use two alternative repre- 
sentations of the baryon-antibaryon amplitude A(pi)A(p 2 ) 
-+ A(p[)A(p' 2 ). 

In the ^-channel representation the amplitude is the 
sum of partial waves in the ^-channel with definite quan- 
tum numbers: spin S, angular momentum L and total 



2.2 Nucleon-antinucleon scattering amplitude 

Let us write the s-channel expansion for a pair of nucleons 
where N = (p, n) forms an isodoublct. The systems pn 
and np have isospin 1=1, and the s-channel expansions 
of their scattering amplitudes are determined by formulae 
which are analogous to those for A A, Eq. (2). The systems 
pp and nn are a superposition of two states, with 1 = 
and 1=1. The nucleon-antinucleon amplitudes read: 

p(pi)n(p 2 ) -> pip'Jnip'z) (1 = 1): 

(pV/i 1/2 , 1/2 1/2) Mi (*,*,«) = Afi (*,*,«) , (4) 
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p(Pi)p(P2) p(p!)p(p 2 ) (I = 0, 1) : 

(Ci 1 /^ 1/2 , 1/2 -1/2) Afi(s,t,u) 

+ (^i°/2 1/2 , 1/2 _ 1/2 ) 2 M ( S ,t, U ) = 

^M 1 (s,t,u) + ^M a (s,t,u) , (5) 
p(Pi)p(P2) -> n(pi)n(p 2 ) (7 = 0, 1) : 

^1/2 1/2 , 1/2 -1/2 Cl/2 -1/2 , 1/2 1/2-^1 ( S >*i U ) 

+ ^1/2 1/2 , 1/2 -1/2 Cl/2 -1/2 , 1/2 l/2 M o(s,i,u) = 

^Mi («,*,«)- ^M (a,t,«) . (6) 

Note that, by writing the NN (or NN) scattering am- 
plitudes, one can use alternatively either isotopic Pauli 
matrices (I/y/2, t/V2) or Clebsch-Gordan coefficients. In 
(4), we use Clebsch-Gordan coefficients which allows us to 
consider reactions in which states with I > 1/2 are pro- 
duced. 

The s-channel operator expansion for NN — > NN can 
be written as 

M I (s,t,u)= J2 (Wi)Q^'iAk'W(-P2)) 

S,L,L',J 
Ml-.-M .7 

x (^(-^Q^W^Pi)) ^'^(s). (7) 

Since the two masses are equal, k = k± holds. In Eq. 
(7), the summation is performed over all states (as well as 
for the A A scattering amplitude). The spin- momentum 
operators Qui. J -uj( k ) f° r tne states with J = 0,1,2 are 
given in Section 4. 



2.3 Amplitude for pA — ► pA scattering 

It is convenient to present the amplitude pA — ► pA pre- 
cisely in the same technique which was used in the con- 
sideration of the s-channel fcrmion-antifermion system. 
To this aim, we declare p being a fermion and A an an- 
tifcrmion. Then, the s-channel expansion for the pA — > pA 
scattering amplitude reads: 

M NA ^ NA { s ,t,u)= ]T (Mp[)Q*tiAk'^ A (-p> 2 )) 

S,L,L' ,J 

x m-P2)Q s ^A k ^( P i)) . (8) 



2.4 Amplitude for AA — > AA scattering 

Let us present the amplitude AA — ► /l/l in the technique 
which was used for reaction pA — > So, we declare 
the 1st yl to be a fermion and the 2nd one to be an an- 
tifcrmion. One can distinguish between them, for exam- 
ple, in the cm. system labeling a particle scattered into 



the backward hemisphere as " antifcrmion" . Then the s- 
channel expansion for the AA — > AA scattering amplitude 
reads: 

M AA ^ AA ( s ,t,u) = £ (Mp[)Q s ^iAk'W A (-P2)) 

S,L,L',J 
M 1 - ■ • M J 

x m-P2)Q*t J .„A k )MPi)) ^Kfto ■ (9) 

In this reaction, a selection rule for quantum numbers 
caused by the Fermi statistics should be taken into ac- 
count, such as: 

= -1. (10) 

Therefore, the following states contribute into (9) only: 

5=1: (7=1;J = 1), (7 = 3; J = 2,3,4), ... 

5 = 0: (7 = 0;J = 0), (7, = 2; J = 2), ... (11) 

2.5 Nucleon-nucleon scattering amplitude 

Nucleon is an isodoublet with components p — > (7 = 
1/2, 7 3 = 1/2) and n -> (7 = 1/2, J 3 = -1/2). The sys- 
tems pp and nn have total isospin 7 = 1, and the s-channel 
expansions of their scattering amplitudes are determined 
by formulae analogous to those for AA, Eq. (9). The sys- 
tem pn is a superposition of two states, with total isospins 
7 = and 7=1. The amplitudes read: 

p(Pi)p(P2) — ► P(p'i)p(p 2 ) (7=1): 

(Ci/2 1/2 , 1/2 1/2) Mi (*,*,«) = Mi(s,t,u) , (12) 
p(pi)n(p 2 ) p(p' ± )n(p' 2 ) (7 — 0, 1) : 

1/2 , 1/2 -1/2 

+ (Cv2 1/2,l/2-l/ 2 ) M>(*,t,«) = 

^Mi(s,t,u) + ^Mo(a,t,«) , (13) 
n(pi)n(p 2 ) -> n(p' 1 )n(p' 2 ) (7 = 1): 

( C 'i 1 /"2 1 -i/2 , 1/2 _ 1/2 ) Mi(s,t,u) = Mi (*,*,«) . (14) 

The s-channel operator expansion gives for Mj(s,t,u) in 
the reaction pn — ► pn (7 = 0): 

M (s,t, U )= 2 (^(Pi)Q^: j „ J (fe , )^(-^)) 

S,L,L>,J 

X (C(-P2)Ql" w (fc)V P (pi)) 4 S ' L ' L ' J \s), 

5=1: (7 = 0; J = 1), (7 = 2; J = 1, 2, 3), ... 

5 = 0: (7 = 1; J = 1), (7 = 3; J = 3), ... (15) 

and for 7 = 1: 

M 1 ( S ,t,u)= (MPl)Q S ^iA k 'Wn(-P2)) 

S,L,L',J 
H 1 ■ • ■ M J 
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5 = 1: (L = 1; J = 0,1,2), (7 = 3; J = 2,3,4), ... 
5 = 0: (L = 0;J = 0), (L = 2;J = 2), ... (16) 

The selection rule for quantum numbers in (15) and (16) 
is caused by the Fermi statistics. 

Analogous partial wave expansions can be written for 
the reactions pp — > pp and nn — > nn (7 = 1), with an ob- 
vious replacing in (16): n — > p for — > pp and p — > n for 
nn — > nn. Here, as for A/l — > A/1, declaring one nucleon as 
a fermion and the second one as antifcrmion, one distin- 
guishes between them in cm. system labeling a particle 
scattered into the backward hemisphere as " antifermion" . 

3 Unitarity conditions and lf-matrix 
representations of baryon-antibaryon and 
baryon-baryon scattering amplitudes 

Here, we write down the unitarity conditions and give the 
A"-matrix representations of the baryon-antibaryon and 
baryon-baryon scattering amplitudes suggesting that in- 
elastic processes are switched off (for example, because 
the energy is not large enough). Generalisation of the K- 
matrix representations in case when inelastic channels are 
switched on can be performed in a standard way. 

3.1 AA scattering 

In this subsection, we consider the unitarity condition for 
the amplitude with J = L. The generalisation for J = 
7±1 amplitude is considered in the last subsection. For the 
amplitude AA — > AA of Eq. (2), the s-channel unitarity 
condition reads for J = L (wc rc-definc A < ; s S ' LL ' J \s) — > 

A aA-^aA^ as follows: 

E Wx)Qlt J .„W)r{-&)) 



X E [(^(Pl)Q&'jW^(-P2)) 



The projection operator 0^h"'^'L is presented in Section 4. 
The phase space is determined as 

2tt 4 

d^iipuPi) = —d<P 2 (pi,P2) = 

-(27Tf 5^>(P-pi -pi) . 3o — , 3o — • 20 
2 (2tt) j 2pio (2?r)' ;i 2p2o 

The projection operator 0^},'"^J, obeys the convolution 
rule, = 2J+ 1 , that gives: 

p [ !a j \ s ) = ^tVt / ^"'M (21) 

xS P (g^^ J (fc")(-p' 2 ' + m A )Q s ^ J (k")(p'; + m A )) . 

The unitarity condition (18) results in the following if- 
matrix representation of the amplitude AA — ► AA 



(22) 



3.2 AA scattering 

Likewise, we consider the unitarity condition for the AA 
scattering amplitude. The s-channel unitarity condition 
for the amplitude AA — > A/1 with J = L reads 

, A (S,LL,J), s _ 1 (SLJ) ( )A (S,LLjy ( )A (S,LL,J) ( , 

1111 a aa-^aa \ s ) — 2 Paa \ s > aa^aa \ s ) a AA^AA \ s )i 

where the identity factor 1/2 is introduced. In this way, 
we keep the definition (20) for d^ip'iiP'i)- 
We have 



(23) 



xS P (Q S »t J .»A k "K-p2 + mA)Q S ^{k"){p'( + m A j) . 
The convolution rule O^ ;^ = 2J + 1, gives us 

P { aa J) {*) = / dS 2 (p'(, p' 2 ') (24) 

xSp {Q S ^A k ")(-P2 + mA)Q S ^A k "M + , 



Finally, one has: 



where 



OH-l-.-Hj (SLJ) 

u ^...^Paa 



(s) = J dMp'lp'i) (19) 



(17) thus leading to identical definitions for p^ AA J \s) and 



P A S A J \ s )i see (21)- The unitarity condition (23) results 
in the following AT-matrix representation of the amplitude 
AA -» AA: 



a(S,LL,J), x 



A' 



(S,LL,J) 
AA^AA 



(s) 



(25) 



Let us emphasize the appearance of the identity factor 1/2 
in the denominator of (25). 
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3.3 Nucleon-antinucleon partial wave amplitude 

The AT-matrix representation for NN scattering ampli- 
tude is written precisely in the same way as for the AA 
case. The only new aspect as compared to AA is that the 
NN scattering is determined by two isotopic amplitudes, 
see (5) and (6), with I = 0, 1: 

pn^pn (1 = 1) : ^Mi(s,i,u) + ^M (s,t,u) , 
pp^nn (J = 0,1) : \m x (s, t, u) - ^M (s, t, u). (26) 

Being expanded over the s-channel operators Qf^„ , (A;)® 

Qf,i.. J iij{k')i these amplitudes are represented through 

partial wave amplitudes A)^' L L,J \s) and A[ S ' L L ' J \s). 
The unitarity condition for these amplitudes leads again 
to a A-matrix representation. 

As above, we consider here the one-channel amplitude, 
first for J = L. The two-channel amplitudes (S = 1, J = 
L ± 1) are presented below in Section 3.4. 

The imaginary part of the amplitude A^ S ' LL ' J \s) with 
I = 0, 1 and J = L satisfying the s-channel unitarity 
condition reads 

J m 4S,LL,J) (s) = ^LL,J) (a)A (S,I,I,,J)- (a)A (5,LL,J) (a)) 

(27) 



where 

(S.LL',.1) 

"nn 



to 



1 

2J+ 1 



d<P 2 (pi,P2) 



(28) 



xS P (0^ L : 7 MJ (fc)(-p 2 + mjvJQ^WGJi + m;v)) . 

The unitarity condition (28) gives us the following A- 
matrix representation: 



A (S,LL,J ){s) = 



K 



(S,LL,J) 



to 



l-ipf^ L ' J \s)Kf^ J \s) 



• (29) 



to those for nucleon-antinucleon scattering, eqs. (27) and 
(28), except for the factor \ in the amplitude originating 
from Fcrmi-Dirac statistics. 



A (S,LL,J ){s) 



to 



K {S,LL,J) ^ 

i^^^F 1 ^ ■ (31) 



(ii) Partial wave amplitudes for S = 1,J = L±l. 

In this case, four partial amplitudes form a 2 x 2 matrix 
given by 



A<F=i,l=j±i,j) {8) 



Ai s^-^j-i,j ){s ^ A (s=ij-^j + i,j)^ 



A (s=i„/+w-i,./ )(sM , 
The A-matrix representation reads 



(S=1,J+1-»J+1,J) 



1 (S=l, i =7±l,7) (s) = ^(5=l,L=J±l,J) (s) 



(^) 



^S=1,L=J±1,J) 

Pnn 



(S=1,L=J±1,J) 



to 



with the following definitions: 



(S=1,L=J±1,J) 



to 



(S=1,J-1-»J-1,J) 
<S=1,J+1->J-1,,7) 



(s), A 
(a), A 



(S=1,J-1-»J+1,J) 
(S=l,J+l^J+l,J) 



to 
to 



^ L = J±1 ' J) (s) 



(S=1,J-1->J-1,J), \ (S=1,J-1->J+1,J)/ \ 

Pjvjv to> Pjvjv ' to 

(S=1,J+W-1„7), s (S=l, J+l— > J+l, J) i \ 

Pnn \ b )> Pnn \ b ) 



Note that the matrices fff 1,L J±1 ' J ) 

J^(S-1,L-J±1,J) ^ gy mme ^ r j Ca J. 



(s) and 



(S=l, J— 1— > J+l, J) 

Pnn 



s ) = Paw 



(S=1,J+1->J-1,J) 



A 



(S=l,J-l->J+l,J)^ s ^ _ ^.(S=l,J+l-»J-l,J) 



(32) 



(33) 



(34) 



to, 



'to- (35) 



3.4 Nucleon-nucleon scattering amplitude 

The pp and nn systems are pure I — 1 states, while the pn 
is a superposition of two states with total isospins 1 = 
and 1=1. The amplitudes read: 



pp^pp, nn — > nn (1=1) : M\(s,t,u), 



(30) 



prwpn (1=0,1) : ^Mi(s,t,«) + ^M (a,t,«). 



Let us emphasize that the definitions of the phase spaces 
for AW and NN systems coincide: Pnn~* L ,J \s) = 
p y N '^ ' ' (s) . In the equation imposing the unitarity con- 
dition (as well as in the AT- matrix representation), the 
identity of particles in the AW systems is taken into ac- 
count directly by the factor 1/2. The unitarity conditions 
for the A A, AA and AW two-channel partial wave ampli- 
tudes for S = 1 and J = L ± 1 are written similarly. 



The expansion over s-channel operators <3^f^ MJ (fc)<8> 
Qfii.- J n /to) ^ s a representation of these amplitudes through 
partial wave amplitudes Aq S ' L L ' J \s) and A[ S ' L L ' J \s). 
(i) Partial wave amplitudes AW — > AW for J = L. 

For J = L, the amplitude A\ S,LL,J) (s) with 7 = 0,1 
satisfying the s-channel unitarity condition are identical 



4 Nucleon-nucleon interaction operators 

In this Section, the proton-proton interaction operators 
are constructed. These operators are constructed using an- 
gular momentum and spin operators, whose properties are 
discussed below. 
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4.1 Angular momentum operators 

The angular-dependent part of the wave function of the 
composite state is described by operators constructed us- 
ing relative momenta of particles and the metric tensor. 
Such operators (we denote them as xj^}..^ , where L is the 
angular momentum) are called angular momentum opera- 
tors; they correspond to irreducible representations of the 
Lorentz group [145,147]. They satisfy the following prop- 
erties [145]: (i) Symmetry with respect to permutation of 
any two indices: 



X (L) 



x (l) 



■ML ' 



(36) 



(ii) Orthogonality to the total momentum of the system, 
P = h+k 2 : 



p x^ L ) — n 



(37) 



The traceless property for the summation over two any 
indices: 



X (L) 



■ Hj-.-HL 



0. 



(38) 



Let us consider a one-loop diagram describing the de- 
cay of a composite system into two spinless particles which 
propagate and then form again a composite system. The 
decay and formation processes are described by angular 
momentum operators. Due to the conservation of quan- 
tum numbers, this amplitude must vanish for initial and 
final states with different spin. The S-wave operator is a 
scalar and can be taken as a unit operator. The P-wave 
operator is a vector. In the dispersion relation approach, it 
is sufficient that the imaginary part of the loop diagram, 
with S and P-wave operators as vertices, is equal to 0. In 
the case of spinless particles this requirement entails 



/ 



dQ (1) 
4tt " 







where the integral is taken over the solid angle of the rela- 
tive momentum. In general, the result of such an integra- 
tion is proportional to the total momentum of the system 
Pft (the only external vector): 



4ir " 



XP, 



(40) 



Convolution this expression with P^ and demanding A = 
0, we obtain the orthogonality condition (37). The orthog- 
onality between the D- and S-waves is provided by the 
traceless condition (38); conditions (37), (38) provide the 
orthogonality for all operators with different angular mo- 
menta. 

The orthogonality condition (37) is automatically ful- 
filled if the operators are constructed from the relative 
momenta fc^ and tensor g-r v . Both of them are orthogonal 
to the total momentum of the system, see eq.(3). In the 
cm. system, where P = (P ,P) = (y/s,0), the vector k 1 - 
is space-like: k 1 - = (0, k). 



The operator for L = is a scalar (for example a unit 
operator), and the operator for L — 1 is a vector, which 
can be constructed from kj^ only. The orbital angular mo- 
mentum operators for L — to 3 are: 



X (0) 



X (2) 
X (3) 

H1H2H3 



4 1} = K 



-I k^k 1 - — - k 2 a 



(41) 



k\ 



k 1 - k^ k^ — 

r (L) 



(9 Hi 1*2^(13 +5a*1M3*>2 +^2M3^Ml) 



The operators Xj Jil '.. tlL for L > 1 can be written in the 
form of a recurrent relation: 



X (L) 



7 a 

L 



k 1 - 7, a 



2L 



L 2 



1(YX^ g ± 

i=l 



i,j=i 

The convolution equality reads 



■ Hj-iUj+i—ULO' 



). (42) 



k 1 - 



k 2 



(43) 



Based on Eq.(43) and taking into account the traceless 

property of xj^}..^ , one can write down the orthogonali- 
ty-normalisation condition for orbital angular operators 

dQ 

\ ' '• ■' in" 
L 



I 



in ^ L (k^)X^l L ,(^) = 5 LL ,a L kf, 



= n 



21 - 1 



1=1 



1 



(44) 



Iterating Eq. (42), one obtains the following expression for 



(39) the operator X^ 



h 1 - h 1 - h 1 - h 1 - h^- 
>1>2>3>4 ' ' • >L 



k\ 



9^,-k±k: 



■ k, 



+9mH3 fc M2 fyi4 • • • ^HL + • • ■ 



+ 



ki 



(2L-l)(2L-3) 
+9, 



( n 1 - n 1 - k 1 - k 1 - k 1 - 



(45) 



4.2 Projection operators and boson propagator 



The projection operator Oj£;;;j?£ is constructed from the 
metric tensors and it has the following properties: 



{1\...{1L V\---Vl J *-Vt...VL 1 

Qfll...fl L QOLt...Ct L = flUl—flL 

KJ a\...a L ^V!...vl ^vx-.-vl " 



(46) 
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Taking into account the definition of the projection oper- 
ators (46) and the properties of the AT-operators (45), we 
obtain 

k^.-.k^O^ =JL X W VL {k ± ). (47) 

This equation presents the basic property of projection 
operator: it projects any operator with L indices onto the 
partial wave operator with angular momentum L. 
For the lowest states, 



= 1, 0$ = 



Q^iPi = I ( g -L Q ±- +q -L _L _^ -L J 

^1^2 2 I ^Pl"! ^P2V2 1 & piV\ t^Hp\P7.av 



.(48) 



For higher states, the operator can be calculated using the 
recurrent expression 



1 / L 

QUI... PL — _[ \ ^ „-L Qpi-.-Pi-lPi+l-.-PL 
^V 1 ...V L L 2 \ ' Vpi V j V\—Vj-\V i + \...Vl. 



(49) 



(2L- l)(2L-3) 



y \ " „-L „-L f)Pl...pi-lp i+ l...p j -lPj + l...pL | 

Z_/ yPiPjUVkfrn ^V 1 ...V k _ 1 V k + 1 ...V m - 1 V m+1 ...V L J' 
i<j / 
A; < ?7i 

The product of two A-operators integrated over solid 
angle (which is equivalent to an integration over internal 
momenta) depends on external momenta and the met- 
ric tensor only. Therefore, it must be proportional to the 
projection operator. After straightforward calculations, we 
obtain 

/ 



% 4^(^)*£U(^=^<« • (50) 



2L + 1 V1 - VL 

Let us introduce the positive value |fc| 2 : 

2 _ 2 _ [s-(mi+m 2 ) 2 ][a-(mi-m 2 ) 2 ] 
1*1 — *±- ^ 



(51) 



In the c.m.s. of the reaction, fe is the momentum of a 
particle. In other systems, we use this definition only in 
the sense of |fe| = \f-k\_\ clearly, |fc| 2 is a relativistically 
invariant positive value. Then, Eq. (50) can be written as 



■ML 
VL ' 



(52) 



The tensor part of numerator of the boson propagator 
is defined by the projection operator. Let us write it as 



4.3 Spin operators of two-fermion systems 

The wave function for fermion particles with the momen- 
tum p is described as Dirac bispinor: 



u(p) 
u(p) 



1 



(po + m)uj 
S^my/po + m V 
(uj*(p + to), -a;* (per)) 
V^m^Pa + m 



(54) 



To construct the operators for the two-fermion system, 
one should also introduce the charge-conjugated bispinors: 



u(-p) 



^y/po + to \ (Po + ™V 



[pa) 



u{—p) = —i 



.(w'*(po-) , -Lj'*{p + m),) 



'2rriyJpo + to 



(55) 



Here, the to and u>' represent 2-dimensional spinors, lo* 
and u'* are the conjugated and transposed spinors. The 
normalisation condition can be written as 

u{p)u{p) = -u(-p)u(-p) = 1 , 
m+p 



u{p)u{p) = 

polarisations 

u(-p)u(-p) = 

polarisations 



2m ' 

-TO + p 

2m 



(56) 



where p = p^"/^ . 

Let us consider a two-fermion system with the total 
momentum P = k\ + fc 2 and relative momentum k = 
(fci — fc 2 )/2, where k\ and fc 2 are their individual momenta, 
P 2 = s. For the sake of generality, let the fermions have 
different masses, mi and m 2 . The two-fermion system can 
form two possible spin state, S = (singlet state) and 
S = 1 (triplet state). The spin operators for these states 
act between bispinor and charge-conjugated bispinor, 
u(— ki)S^u{k2) and have the following form: 



where 



»75 



7 U 



(57) 



(58) 



It should be noted that u(— ki) and u(fc 2 ) have opposite 
parities, so u(— fci)75«(fc 2 ) is a scalar and u(— fci)u(A: 2 ) is 
a pseudoscalar. 

As is shown below, that the 7 M operator leads to the 
mixture of states with total momentum L + 1 and L — 1. 
So, let us introduce the operator for the pure 5 = 1 state: 



s& =r ± = 



^kj 



M{yfs + M){s-5 2 ) 



(59) 



JTPl—PL _ (_-\\L QP1..-PL 
^ Vi. ..VL V Vi ...VL ' 



(53) 



This definition guarantees that the width of a resonance 
(calculated using the decay vertices) has a positive value. 



where M = mi + to 2 and 5 = mi — m 2 . In the nonrela- 
tivistic limit, this operator is equal to the spin-1 operator 
<r and satisfies the orthogonality of the triplet states with 
the same parity. 
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4.4 Operators for 1 L J states 

In case of a singlet spin state, the total angular momentum 
J is equal to the orbital angular momentum L between the 
two particles. The ground state of such a system is 1 S'o 
( 2S+1 Lj) and corresponding operator is just equal to the 
spin-0 operator of Eq.(57). For states with orbital 
momentum L, the operator is constructed as a product 
of the spin-0 operator and the angular momentum 
operator X^...^,: 



V^^^-^i^X^..^). (60) 

The normalisation factor which is introduced here simpli- 
fies the expression for the loop diagram (see below). 



4.5 Operators for 3 Lj states with J — L 

The ground state in this series is 3 Pi, so one should make 
a convolution of two vectors, Sjx^ and x'jp that creates 
a J = 1 state (vector state). In this case the vertex oper- 
ator is equal to e VlV f 1 ^ ri k^P 1 . For states with higher or- 
bital momenta, one needs to replace k^~ by xi' 1 ^ and 
perform a full symmetrisation over v\ , v 2 , ■ ■ ■ , vj indices, 
which can be done by a convolution with projection oper- 
ator O^ytfL . The general form of such a vertex is hence 
given by 



V L=J ~f c 'v xi J ^ p nm—t*j 



(61) 



Using eqs.(45) and (47), one has 



2J- 1 



(62) 



Finally, using Eq. (167) the vertex operator can be written 
as: 



V, 



L = J 



(2 J + 1) J ie avil - 1v kfP 1 Z° 



(J+l)aj 



(63) 



where normalisation parameters are again introduced. No- 
te that due to the property of the antisymmetric tensor 
e ar jfj the vertex given by Eq. (63) does not change, if one 
replaces j n by the pure spin operator r v . 



4.6 Operators for 3 Lj states with L<J and L>J 

To construct operators for 3 Lj states, one should multiply 
the spin operator j a by the orbital momentum operator 
for L = J + 1. So one has: 



V L<J ~ ->„ A (J_1) 



Using Eq.(167) from Appendix 3, we write the vertex op- 
erator in the form 



and for the pure spin operator as 



■TrL<J _ p j " 



(65) 



(66) 



The normalisation constant is chosen to facilitate the cal- 
culation of loop diagrams containing such a vertex. 

To construct such an operator for L > J one should 
reduce the number of indices in the orbital operator by a 
convolution with the spin operator: 



and for pure spin state: 



\rL>J p y~ 



J+l 

aj 



J+l 

aj 



(67) 



(68) 



5 Calculation of the NN — > NN amplitude 
5.1 Structure of the amplitude 

Let us consider N(qi)N(q2) — > N(ki)N(k 2 ) transition am- 
plitude with q 1 , q 2 , k\, k 2 being the nucleon momenta, and 
k = (ki — fc 2 )/2, q = (qi — q 2 )/2. In this section the struc- 
ture of such amplitude for different initial and final states 
is derived. 

We start by considering NN — > NN amplitudes for a 
singlet state. For the 1 Lj state with L — J, the amplitude 
is given by 

A s = u{-q 2 )V^... IJlJ {q) u{q x ) 

x F£"ff u(ki)V Vl ... K ,(k) u{-k 2 ) . (69) 



For sake of simplicity, we omit here and below the invari- 
ant part of the amplitude, which will be considered later 
on. Using Eq. (60), the amplitude reads 

A s = -u(-q 2 ) 75 u(qx) u(ki) 75 u(-k 2 ) 

x(\k\\q\r(2J+l)Pj(z), (70) 

where z = (kq) / (\k\\q\) is the cosine of scattering angle 
in cm. system. 

The transition amplitude for the triplet state 3 Lj with 
L < J is the following: 



A L<J 



H~q2)v^%,(q) u( qi ) 



V1...UJ 



(64) 



x Fg;;;% u{kx)V^< J VJ {k) u(-k 2 ) . (71) 
Using Eq. (167), this amplitude can be written in the form 

Af< J = f^ fi<H{\h\\<l\) J -\ (72) 
i=i 
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where 



h 


= u(- 


-qz) 7m u (h) u ( k i) 7m u(-k 2 ) , 


h 


= «(- 


-q 2 ) q u(qi) u(ki) q u(-k 2 ) , 


h 


= «(- 


-q 2 ) k u(qi) u(fci) k u(-k 2 ) , 


h 


= «(- 


-q 2 ) q u(qi) u(ki) k u(-k 2 ) , 


h 


= «(- 


-q 2 ) k u(qi) u(fci) q u(-k 2 ) . 



(73) 



Here 
ai = 

a 4 = 



P'j(z) 



a 2 



J|q| 2 



«3 



1 



(P;_ 2 (z)-2P' 7 _ 1 (z)) , a 5 



P'J-i(z) 
J|fc| 2 
P'J{z) 



.(74) 



J|fc||<z| * J — ' — " ' J|fe||g| 

Likewise, the transition amplitude for triplet state 3 Lj 
with L > J is as follows: 



j+i 



(75) 



where 



«i = -- 



P'j+i(z) 



i=i 



a 2 = a 3 = 



j + i 



a4 



J + 



J + i 

1 - (P'J(z) + (2J+l)P' J+1 (z)) , a 5 



j+r 

If the spin-1 operator is defined as *y v , there is a mixture 
between two triplet amplitudes with L > J and L < J. 
The corresponding transition amplitudes are given by 

V^%Aq) (-1)' 7 0%;;% V v L > J VJ {k) = 



Ifel^kl- 7 - 1 , 



^» (-1)' 7 0%;;% VZ< J VJ (k) 



E /» a, |fc| J ^ql 



(77) 



i=l 



where 
ai = - 

«3 = 



J PM 

J+12J-1 



J p'UM 

J+l 2J-1 



«2 = 



&4 = as 



J ^-l(^) 

J+l 2 J - 1 



(78) 



J+12J-1 

If one uses the operators based on the pure spin-1 operator 
given by eqs. (66) and (68), the functions fi, f 2 , ■ ■ ■ fs are 
substituted by the new functions /i, f 2 , . . . f 5 as follows: 

li = UMn , (79) 

where the transition matrix Mji is equal to 

/ 1 \ 

K -jy- -/^/tV) 



— K 







VI n 

M U 

y/sK.(k q ) S 



V 



M 

o 



M 2 





-/t(/rV) 
2 (fc^) 2 



(80) 



Then, the transition amplitudes for the 3 Lj triplet state 
with L < J and L < J are 



^ <J = E /^a, (IfelM)- 7 " 1 , 

Af> J = E /iM^Oi (|fc|| 9 |) J+1 • (81) 

The transition amplitude for the 3 Lj triplet state with 
L = J is given by 

Af=^S(^ 2 )^ 7 (a) U ( gi ) 

x F%;;;% uih^i^k) u(-k 2 ) . (82) 

Using expressions given in Appendix 3, this amplitude can 
be written in the form 

A L = J = (f l0l + f 5 a 5 + f e a e )(\k\\q\) J , (83) 

where 



(76) / 3 = u(-q 2 ) 7 P u( qi ) u(fci) j u u(-k 2 ) n v , (84) 



2J+ 1 



a 5 



(J + l)J|fc| |q| 
5.2 One-loop diagrams 



(J+l) J 



The calculation of one-loop diagrams for different vertex 
operators is an important step in the construction of a 
unitary NN amplitude. Let us start from the loop dia- 
gram for the singlet state and derive all expressions for 
the case of different particle masses, m\ and m 2 . Taking 
into account that 

75(^1+^1)75(^2-^2) =2(s-(5 2 ), (85) 

where 5 = mi — m 2 , the one- loop diagram for the singlet 
state is equal to 



-/ 



df2 

47T 



Sp 



V^ 1 ... flJ (k ± )(m 1 + k)V Vl ... VJ {k ± ){m 2 - k 2 ) 



2(s-S 2 )\k\ 2J 0^ j '(-l) J - 



(86) 



The factor (—1) is related to the fermionic nature of the 
baryon in the loop. 

To calculate one-loop diagrams for different triplet sta- 
tes, the following relations are helpful: 



Sp 7 M L (mi + fci)7^(™ 2 -fc 2 ) = 2{s-5 2 )g^+Sk^ 



Sp 



r„(mi + fci)r„(m 2 -feO =2(s-o 2 ) 3 ^. 



(87) 



Using these relations and eqs. (172)-(178) given in Ap- 
pendix 3, we obtain the following results for the L < J 
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and L > J states: 
dQ 



I 



4tt 



Sp 



V^iAmi + hW^imz - k 2 ) 



I 



2(s -^ 2) -^Tl) l fc l 9(J-1) °«/(-i) J . 

fspf^^K + fcO^^ma - fc 2 ) 



2(s-<5 2 



8(J+l)|fe| 
2J+1 



|fc| 2(J+1) OK/(-l) J (88) 



In case of spin-1 operators, the two triplet states with the 
same parity are not orthogonal to each other; the inter- 
ference loop diagram is equal to 



df2 
Air 



Sp 



V^iMi + ki)V v \> J VJ {m 2 - k 2 ) 



L>J 



8^^|fc| 2(J+1) 0^/(-l)^ 



(89) 



The one-loop diagram for the L = J triplet state is given 

by 



-7-Sp y„ 1 ... w (ffli + ki)(k^)V Ul ... Vj (k^)(m 2 - k 2 ) 

'ITT 



= 2{ s -8 2 )\k\ 2J o^:::^(-i) J . 



(90) 



Direct calculations show that the transition loop diagrams 
between the triplet state with L = J and the triplet states 
with L > J and L < J vanish. For vertex operators de- 
scribing pure spin states (66) and (68), one has the fol- 
lowing one-loop diagrams: 



-/ 



df2 

47T 



Sp 



2{s-s 2 )\k\ 2{J+1) o<ii:X : ;{-i) J . 



(91) 



5.3 Cross sections 



The expressions for one-loop diagrams can be used to cal- 
culate cross sections for different spin-orbital momentum 
states. The cross section is given by 

Mq\Vs \q\ 16tts J 4tt 1 1 v ' 



\q\ 16tts 



To calculate the amplitude squared \A\ 2 , one can use the 
expressions for the one-loop diagram given by eqs. (86), 
(91), (89) and (90). 

For the 1 Lj state, one has 



da = 



2J+1 
647rm 2 



| fc |2J+l 1^2.7-1 s 



(93) 



For the 3 Lj state (L = J), the result is 



s . 



(94) 



The decay of the 3 Lj states with [L < J) and 3 Lj (L > J) 
is determined by the sum of two vertices: 



the cross section is equal to 

da = Xfdau + X\da 22 + Xi\ 2 (dai 2 + da 2 i) 



(95) 
(96) 



where 
dan 



2J+ 1 
2567rsm ; 



2 M 2 '- 1 \i\ 2J - 3 



X 


L_sm 2 - 








2J+1_ 




2J+1 



d(J22 = 2 _ J+ \ \ k \ 2j + 3 i-i 2j+i 



2567rsm 2 



2s 



8(J+l)|fc| 



2J+ 1 



2s- 



8(J+l)l9| S 
2J+ 1 



da- 



21 



_i J< y J + jj |q| 2J + 3 |fc| 2J+1 



47rsm 2 2 J + 1 



(97) 



For pure spin-1 operators, V^< J and V^ > J the cross 
section reads 



where 



dan 
da 22 



da = \\da 11 + \\da 22 , 

_ (2J+ l)s | fc |2.7-l |g|2J-3 

647rm 2 

_ (2J+ l)s ,2,7+3 |„|2,/+1 

_ 64W 191 



(98) 



(99) 



6 Decay into 3/2+ and 1/2+ particles 

Let k\ be the momentum of the 3/2+ particle and k 2 the 
momentum of l/2 + . In this case, there are two spin states, 
5=1 and S = 2. Let us start from the S = 1 states. Such 
states are constructed using the vector spinors ip a for 3/2- 
spin particle and spin operators. 

For 3 Lj (J = L — 1), that corresponds to 
(0~, 1 + , 2~, 3 + . . .) states, and the operators read 



(100) 



For 3 Lj ( J = L + 1), the operators are given by 

vW%=ii 5 xii:Zo%;;% , (101) 
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where the projection operator is needed for index sym- 
metrisation. For 3 Lj (J = L) (1", 2+, 3", 4+ . . .) , the 
operators can be expressed as 

= ^e aiKv hP v X^-ZO^ . (102) 

In case of S = 2, there are five operators. For 5 Lj 
( J = L + 2) the operators are given by 

wgln = ^ 1 (fci)7a 2 o^ 2 2 ^3 J :^^i::-;>(- fc 2). 

= l a2 O^X^,0^;;Z , (103) 

for 5 Lj ( J = L - 2) by 

= l a .Xi{tf Vl ... VJ , (104) 
for 5 Lj (J = L) by 

V£*„ = l P O^X^ Mj O^;;Z , (105) 
for 5 Lj (J = L - 1) by 

= 

^ 1 /3r J) fc r p^o«f^ ai (fc 1 )7a 2 ^S... I/ ,op i ;:;>(-fc 2 ), 
^A^^^^^^r 2 ^^^.,.^^;;;^ , (ioe) 

and for 5 Lj (J = L + 1) by 

^(8) = 

=i^/JrA^O^ aa 'ya a ^ J :^0^ 1 ;;7 J . (107) 

The one-loop diagram amplitudes for the correspond- 
ing operators are calculated in Appendix 4. 



7 Example: amplitude for the reaction 

PP — *■ pK + A 

Let us start from scattering with the production of a 
resonance i? in the intermediate state which decays into 
K + A. The diagram for the process is shown in Fig. 1. 
Consider the partial wave amplitude for the pp having 
quantum numbers J = n, L and S in the initial state. 
The general form of the angular dependent part of this 
partial amplitude is 

(u(-fci)Q£'^>(fc 2 )) u(q 3 )N ai ... am (R^KA) 

*^::Z(Q* + to)Qtt? m ^... Vn <-*) 

-{hi fc 2 } , (108) 




y K + (q 2 ) 



A(q 3 ) 



Fig. 1. Reaction pp — > pK + A: pp scattering with production 
of a resonance 7? in the intermediate state 

where P = q\ + q 2 + q 3 = k\ + k 2 . The resonance R with 
spin J = m + 1/2 is produced in the intermediate state 
and decays into a final state meson and a nucleon. 

The initial pp state operator Q^f.'^P is defined by Eq. 
(60) for S = and eqs. (63), (65), (67) for S = 1. If the res- 
onance in the intermediate state has the spin 1/2 (m = 0), 
the same expressions define the Rp state operator. For the 
spin-3/2 resonance in the intermediate state, the operator 

Q ( 0^.0^... Vn is defined by eqs. (100)-(107). The opera- 
tors for the R — > 0~ + 1/2+ transitions were defined in 
[147]: 



N + = X (n) 



1/2", 3/2+, 5/2"... 
N-...^ = il v ^Xi n +% n 1/2+, 3/2-, 5/2+ . . . (109) 

Let us write down the amplitude for the 1/2+ reso- 
nance in the intermediate state. In this case, one finds 

S,L,J 
. <?2 + 93 



X ■ 



S ^BW(s 23 )QifXM-Q^ L ' J Hs,s 23 ) 

(110) 



<?23 



- {h ^ k 2 } , 

/ \ / Q23 fJ ,Q23i>\ „ 
= (52 - q3)v[g»v , Q23 = 92 + 

v So.3 ' 



93 



where A^ s - L ' J \s, s 23 ) is the partial amplitude forpp — > Rp 
and BW(s 2 3 parameterises the resonance R). 

Another type of processes, which may contribute to 
this reaction, is the t-channel exchange of pseudoscalar 
and vector particles (ir and p). This diagram is shown in 
Fig. 2. 

First, consider the exchange of a p meson. The vertex 
operators for the transition baryon — > vector meson + 
baryon are given in [147]. Thus the pp — > R operators are 
given by 



ALL = u{Q23)V {l -^{k^)u{k 2 ) Ptl 



lower 
K 2n 

where 



i(fc 2 - k t ) u {g^ 



Q23nQ23v " 
S 2 3 



2 = 1,2 , 
(111) 
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Fig. 2. Reaction pp — > pK + A: t-channel exchange diagram. 



(112) 



Here, k t = Q23 — k 2 is the p-meson momentum. In case of 
a spin-1/2 resonance in the intermediate state, one should 
use Eq.(112) with n = 0. For the upper operator, one has: 



= u(qi)V^(qtHkl) Pt , , 

5i„ = 2 («i - fc *)^ \9i>» ) 

Summing over the polarisations yields: 



^2 P a PP = ^9ap + 

polarisations 



hah/3 



i = 1,2 

(113) 
(114) 



Finally, we arrive at the following amplitude for the p 
exchange: 

Af = uiq^V^iqiMhMq^N-iq^) 



<?2 + <73 + \/S23 



BW(s 23 )V^(k^)u(k 2 ) 



(-^ + ^), i = l,2. 



(115) 



In case of t-channel exchange of a pseudoscalar meson, it, 
one should substitute the operator V^ l ~^ by N~ , so we 
have 

A„ = u{q 1 )N-{qi)u{k l )u{q i )N-{qi s ) 



g2 + <?3-hv/£23 



BW(s2 3 )N-(k 2 L )u(k 2 ) . (116) 



8 Triangle-diagram amplitude with 
pion-nucleon rescattering: logarithmic 
singularity 

In the amplitudes describing production of three-particle 
final states, the unitarity condition is fulfilled automati- 
cally when final-state rescattering is properly taken into 
account. However, rescattering may lead to singularities 
where the amplitude tends to infinity. The triangle dia- 
gram with the A in the intermediate state gives us an 
example of this type of the process: it has logarithmic sin- 
gularity which under certain conditions (^/s ~ + m A ) 
can be near the physical region. 



Because of yfs <~ mpf + m A , we consider the amplitude 
pp — > NA with L' = (the produced NA system is in the 
5- wave) . The quantum numbers of the final state are then 
restricted to 



J p = l+,2+. 

The initial pp system (1=1) has 

5 = 0: L = 0,2,4,... J p = 0+ 2+ 4+ ... 
5 = 1: L = 1,3,5,... J p = 0-,l-,2-,3-, 

We thus consider the transition 

pp (5=0, L=2, J p =2+) -» NA (5'=2, L'=0, J p =2+) 
The corresponding pole amplitude reads: 



(117) 



(118) 



tpole 
NN—>NNtt 



^NN^NN-n^pp^NA 



/ ~,(S=0,S'=2,L=2,L'=0,J=2) , \ 



V ^ - - im A r A 

x (u(- P2 )i 15 xii](k)u( P1 )y 



lu'U{-p' 2 ) 



(119) 



Here, the factor C^n^nn^ refers to the isotopic Clebsch- 
Gordan coefficients, and 



K lp — 9pp' Pip,' 



(120) 



The numerator of the 3/2— spin fermion propagator is 
written in the form used in [146,147]: 



2M A 

9pvlv 



9^ 



(121) 



The decay vertex g A is determined by the imaginary part 
of the loop diagram A — > Ntt — > A. For the sake of sim- 
plicity, we change in (119): r v >(k' ± ) — > 7^; however, using 
definition (59) one can easily rewrite Eq. (119) in a more 
rigid form. 

Taking into account the rescattering process in the am- 
plitude (119), irN — > Zi — > 7rA r , one has the following 
triangle-diagram amplitude (see Fig. 3): 



a triangle ^triangle 

NN — >NN7r — '-'NN^NNTr'-'pp^NA 



n (S=0,S'=2,L=2,L'=0,J=2), n 



( a «»[/|S 



(2tt) 4 - fc2 - »0 



-f>2 



m zi — Pzt — irn A l a m ~ P 2 



iO 



xg A k 2 1 - - ■ - • — 

^ m\ - p A - im A r A 



9Ak' 2 ^]u{-p' 2 )) 



(u(-p 2 )t 75 X%,(k)u( Pl )) . 



(122) 



Here, 



K lp' ~ 9^'a Plot ) 



>''-LP4 
v 2p" 



9n" a P2a ) 



(123) 
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and 



PA = Pi + fc 7Tl VA = V'l +PTT = P2 + fc 7T , 

P=Pa+P2- (124) 



p(p 2 ) 

p(p" 2 ) A(p * } 




p(p' 2 ) 



Fig. 3. Triangle diagram with final state pion-nucleon rescat- 
tering 



One can simplify (122) by extracting the numerator in 
the singular point that corresponds to 

m 2 A =p'l, m 2 =pf, ml = kl. (125) 

Then, Eq. (122) reads: 

^triangle ^triangle ^(S—O.S —2.L—2.L — 0,J— 2)/ \ 

A NN-+NNtt — NN^NN-it pp^NA \ s ) 

x (s(p' 1 )^fc^(tr)Z\ Al v(p' zl (tr))7^(-p' 2 '(tr) + m)g A 

x (u(- P 2)i~/5Xi 2 J,(k)u( Pl )) 

d A K 1 1 

i(27r) 4 ml~kl- iO m 2 - (p A - K) 2 - iO 

' (126) 



(P - pa + K) 2 - imA^A 



The momenta fc^f^(tr), p^(tr), p£(tr), k'£fe(ti) obey 
the constraints (125). The integral in (126) corresponds to 
the triangle diagram with spinless particles. Its calculation 
is performed in Appendix 5. 



9 Box-diagram singularities in the reaction 

NN -> AA -> NNtttt 

The primary aim of a partial wave analysis is to extract 
the pole singularities of amplitudes, thus determining res- 
onances. Of course, the existence of other singularities like 
threshold singularities should be taken into account. This 
is possible using the if-matrix technique, see [150-152] 
and references therein. Singularities due to resonances in 
the intermediate state need more sophisticated treatment. 

The existence of triangle-diagram singularities, which 
may be located near the physical region of a three-particle 



production reaction, was proven in [153,154]: these singu- 
larities diverge as ln(s — so). Stronger singularities (with a 
(s — so) -1 / 2 behaviour) are related to box diagrams [155, 
156]. 

Here, we present box-diagram and triangle-diagram 
singular amplitudes for the reaction AW — > AA — > NNtttt 
taking into account the spin structure in a way which al- 
lows us to include these singular amplitudes into partial 
wave analyses (this was not yet done in [155,156]). 

Let us introduce the following notations for the two- 
pole and box diagrams in the reactions NN — > AA — > 
NNtttt (see Figs. 4 and 5). 

The initial state momenta are: 

P 1 +P 2 = P, P 2 = W 2 , \{Pi-P 2 ) = q (127) 
Final state momenta: 

(Pi +P3) 2 = S13, (pi +P3 +P2) 2 = s 4 , 



Pl+P3 = fel, 
(P2 +P4f 



■Pz) 



J-fci 



_L/Cl 

Pi 



S24, 



(P2 



P2 H 
(Pi 



Pi 



Pi) 



k 2 



\(P2 



PA+Pl) = Si 

±fe 2 



■Pi)' 



±k 2 
P2 



Xfei 
-P 3 > 



±fe 2 

-Pa 



Pl+P2=P- 



(128) 



Here, the symbol J-ki means the component of a vector 
perpendicular to kf. 



Fp 



Pv 



k 



{hp) 



(129) 



9.1 (NN)s- 
diagram 



state with J p — 0+, two-pole 



In pp collision with 1=1, the S*-wave A A state is pro- 
duced. First, consider the two-pole diagram of Fig. 4. The 



\*( Pl ) 




P(P 2 ) 



Fig. 4. Pole diagram for reaction NN — » AA — » NNtttt 



amplitude for the production and decay of two Z\-isobars, 
NN — > AA — > NNtttt, omitting charge indices and corre- 
sponding Clcbsch-Gordan cocfhcicnts, reads: 



l NN^AA^(Nir)( L Nir) 



(u(-P 2 )u(P 1 ))g nn ^aa(W) 
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x (u(P3)9Ak 



^ M 2 A - s 13 - iM A r A 
A v , v (-k 2 ) 



Ml - S24 - iM A r A 



{-)ki v g A u{- Pi )). (130) 



9.2 Box-diagram amplitude with pion-pion rescattering 

The box-diagram amplitude with pion-pion rescattering 
in the Feynman technique (see Fig. 5) is equal to 



P(P) 




, ji(k ) , , 



£ 7T(k ) X..^ 



Fig. 5. Box-diagram with pion-pion rescattering 



L NN^AA^NN+(wn— ►7T7r) s 



iS-wave/ \ 



x G 



NN^AA (W)(u(-P 2 )u(P 1 ))(u(p 3 ) 



d 4 k' 

i(2TT) 4 



g A k%A tlv ,(k' x )A v , v (~k 2 )(-)k' 2 \,g A 
(Ml - si 3 - iM A r A )(Ml - s' 24 - iM A r A ) 
1 



(ml-k^-iOXml-k^-iO) 



«(-P4))- (131) 



The factor Af"^™ e (s) is the S- wave 7T7r-scattering ampli- 
tude. Here we take into account the low-energy 7T7t inter- 
action only. In the if-matrix representation, it is written 
in the form 



\ S~ wave 

7T7T — >7T7T 



(s) = 



K(s) 



1 / s — 4m 2 

l-ip{s)K{ S y p(s)= i6^v 7^ 



(132) 



In (132), we take into account the full 5-wave as observed 
experimentally including the so-called sigma-meson, inde- 
pendently of its existence. Generally speaking, it is possi- 
ble to account for higher waves as well, but the box dia- 
gram with two A's leads to singularities near the physical 
region of the production process at \fs < 0.6 GeV only. 

The approximation used in the calculation of the box 
diagram (131) is related to the extraction of the leading 
terms of the singular amplitude. To this aim, we fix the 
numerator of the integrand in the propagator poles by 
setting 

kf Ml , k'Z Ml , kl ml, k 2 2n ml, (133) 



which leads in (131) to the substitution 

K ~* fc lU b0X ) = "PS MbOX) , K fc^box), 

k' 2 i fc 2 i(box) = -p 4 ±fc2(box) , k' 2 fc 2 (box)). (134) 
Now, in the cm. system, the momenta fc a (box) read 

fei(box) = (W/2,0,0, ^W 2 /A-M 2 A ), 

fc 2 (box) = (W/2,0,0,~yJw 2 /4~ Ml). (135) 

Here, we denote the four-momentum as k — (ko,k x ,k y ,k z ). 
Under the constraints of Eq. (133), the numerator of the 
integrand does not depend on the integration variables, 
and it can be written separately for the leading singular 
(LS) term: 

A NN-*AA^NN+(nn^nn) s = A ^Z7n ( s )Gn N->Aa(W) 

x (u(-P 2 )u(P 1 )) (u(p 3 ).94(-^ l(bOX) )^'(fci(box)) 

xZ\ 1/ /, y (-fc 2 (box))p 4 L j £(box) g zl u(-p4)^ 
d 4 k' 1 



J if> 



) 4 (Ml - (\ P + k' + P3 ) 2 - iM A r A ) 
1 

x (Ml - (\ V -k' + p A ) 2 - iM A r A ) 

1 

X (ml - (\p + k') 2 - i0)(ml - (\p - k') 2 - i0) 
where 



, (136) 



^p + k' = k l7r , ^p-k' = k 2v , Pi+P2=P- (137) 

The box-diagram integral of Eq. (136) is calculated in Ap- 
pendix 6: in this Appendix, we demonstrate the effects of 
the box diagram on the irn spectra. 



9.3 Box-diagram amplitude with pion-nucleon 
rescattering 

In the Feynman technique, the box-diagram amplitude 
with pion-nucleon rescattering in the resonance state (/ = 
3/2, J = 3/2) reads (see Fig. 6): 



it(P,)/' 




••• n(k' J . *- 
\ v tc(P,) 

P(P 2 > P(P 4 )' 

Fig. 6. Box-diagram with pion-nucleon rescattering 
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AnN^AA^Ntt+(Ntt^Ntt) a — GnN^Aa(W) 

x (u(-P 2 )u(Pi)) 

d ' k ' \{^-K„)r^ J iN _t mH 



i(27r) 4 2 



'AT 



X (M 2 - £f - iM A r A ){M 2 A - k' 2 2 - iM A r A ) 

X Jm^~h~ <0)]^" ( ~ P4) ] ' (138) 
where p^ = p 2 + P3- By fixing the numerator of (138) at 

M 2 A , k' 2 2 - Ml , fc 2 ^ ^ m 2 j ^ ^ m 2^ (139) 



_Lfci 



we write the leading singular (LS) terms of the box-dia- 
gram amplitude as follows: 



A 



(LS) 

N N — *AA — *Ntt-\-(Ntt — >Ntt)a 



G 



NN^AA 



(W) 



X (u(-i^)u(Pi)) 

x^fc^box) -p x - k 2 (box) + P4 )^F A 
xg A (fci(box) -pi+ gA pf^ 1 (b ° x) 
xZ\ MV ,(fc 1 (box))Z\ !y , y (-fc 2 (box))pi iy fc2(box) «(-p 4 )) 
d 4 ^ 1 



i(27T) 4 (m 2 , - (p 4 - fc 2 ) 2 _ z0 ) 



(M 2 - (p 4 - K +Pi) 2 - iM A r A ) 



{Ml - {K + Pa) 2 - iM A r A ){ml - fc 2 - iO) 



• (140) 



The box-diagram amplitude with the pion-pion rescatter- 
ing is given by 

A NN->AA^NN+(™^™)s = A ™ZZ e ( S ) G NN^AA(W) 

x (u{-P 2 )X^Uq)u{Pi)) 
x(u(p 3 )g A (-p^ l{hoK) )A^(ki(box.)) 

x A v « v ( - k 2 (box) )p4 I/ fe< b ° x) g A u ( -p 4 ) ) 

d 4 fc' 1 

z(2tt) 4 (M 2 - (|p + fc' + p 3 ) 2 - iM A r A ) 
1 



(M 2 - (ip - fc' + Pi y - iM A r A ) 
1 

(to 2 - (±p + fc') 2 - i0)(m 2 - (±p - fc') 2 - iO) 



• (142) 



In the leading singular-term approach, the box-diagram 
amplitude with the pion-nucleon rescattering can be writ- 
ten in the form 

A nn-^aa^N7t+(N7t^N7t) a = G N n^aa(W) 

x(,(p3 )g 4(p 2 -P3)^ Mi _^l^ 

x^(fci(box) - pi - fc 2 (box) +p 4 )^ Pzi ffzi 
x (fci(box) - pi +mAr^ ^p^ [t fe , l(box) Zi ([1 / I/ /(fci(box)) 
xZ\ 1/ // l ,(-fc 2 (box))p4j C2(box) g A u(-p 4 fj 
d 4 K 1 



i(27r) 4 (m% - ( PA - fc 2 ) 2 - iO) 



(M 2 - (p 4 - fc OT + pi) 2 - »M 4 r 4 ) 



(M 2 - (fc* + p 4 ) 2 - iM A r A ){ml -kl- iO) 



• (143) 



9.4 (NN) D _ wave state with J p = 2+, two-pole and 
box diagrams 

The production of Z\Z\ near the threshold in the S-wave 
leads to a J p = 2 + state as well and, correspondingly, to a 
strong box-diagram singularity in this wave. In the J p = 
2+ wave, the transition (NN)n-wave — * {AA)s- W ave is 
related to the two-pole amplitude 

A(NN) D ^(AA) s ->(Nir)(Nir) = GnN^Aa(W) 

x(u{-P 2 )X^ v „{Q)<Pi)) 
( — / \ ,± A.^{ki) 

^-fX, '-^ ' <141) 



10 Conclusion 

We have developed a new method for the partial wave 
analysis of data on the baryon-baryon and baryon-anti- 
baryon collision. The method is based on the operator de- 
composition approach which was successfully applied be- 
fore to a number of meson-induced reactions. The article 
emphasises the analysis of reactions with three or four 
particles in the final state, where triangle and box singu- 
larities might play an important role. A full set of partial 
wave amplitudes is constructed for nucleon-nucleon elas- 
tic scattering and for NA and AA production. With these 
amplitudes, expressions for partial widths and for reaction 
cross sections are presented. Some examples how to cal- 
culate contributions from triangle and box diagrams in 
simple cases are explicitly given. The application of the 
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methods developed here to the analysis of new data ob- 
tained and expected from COSY should provide valuable 
information about the hadron spectrum and properties of 
hadron interaction. 
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Appendix 1 

The baryon wave functions ip{p) and ip(p) = ip + (p)jo obey 
the Dirac equation 

(p - m)ip(p) = , ${p){p - m) = . (144) 

The 7-matrices were used in the form 



7o 



/ 
-/ 



7 



a 
-<x 



75 = «7o7i7273 = 



/ 
/ 



To" = 7o , 7 + = -7 , 
with the standard Pauli matrices 



(145) 



1 



-i 



1 



a> 1 1 j • i / o )> a3 v." i, 

a a a h = iSabc^c • (146) 
The Dirac equation gives four wave functions 



, , ifj 

j = 1,2: ipj(p) = VPo + m jap) 



po+m ^3 

• \;(p) \ IK) + hi ( <pf,-<p 
((TP) 



(147) 



+ ■-+ ( g P) 



7=3,4: ipj(—p) = Wpo + m ( po+m Xj 

Xj 



(148) 



ipj(-p) = -WPo + m Xj 



Po + m J 



where <pj and Xj are two-component spinors 



(149) 



which are normalized as follows: 

<Pt<Pt = 6 jt> xtxe = $]£ ■ 



(150) 



The solutions with j = 3, 4 refer to antibaryons. The cor- 
responding wave function is given by 



j = 3,4: r 3 (p) = Cjj(-p), 



where 



C = 7270 



II -a 2 

-cr 2 



The equation (151) reads: 
i = 3,4: 

-cr 2 



^(p) = ( _ CT2 J WPo + m 



( i 
-i 

\-i 0, 



(o- T P) v * 

Po+m X-j 



(151) 



(152) 



(153) 



~WPo + rh (crp 



<JiXj 



Po+m 



<T2Xi 



VPo + m (p-p) , 



po+m rj 



In (153), we have used the commutator 

er 2 (<X T p) = —&1P102 — (JlPlUl — (JzPzCf2 = —{<Tp)(J2- 



(154) 



We define the two-component spinor for antibaryons as 
'0 -1 



<Pj = -*°2Xj = I ! o I Xj 



~Xj2 

X*jl 



(155) 



The wave functions defined in eqs. (147)-(148) are nor- 
malized as follows: 

j,£ = 1, 2 : (tpjip^eXp)) = 2m 5 je , 

jj = 3,4: (^(~ P )M-P)) = -2m 6 jt . (156) 

They obey the completeness relation 

j=l,2 

X] ' t Pja(-p)i J 3fj(-P) = {p- m )a[3- (157) 
.7=3,4 



Appendix 2 

(i) The S'-wave terms in the the nonrelativistic 
limit. 

We consider the operators with L = from eqs. (60) and 
(66) in the cm. system (p 1 = — p 2 = k and p\ = — p' 2 = 
k'). For L = 0, we have the following operators in the 
nonrelativistic approach: 



Q oou (A0 = *T5 



Q 101 (k) = r^ 



I 

1 

er 

-a 



(158) 
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In the cm. system, we have: 



j,f = 1,2 : ip Nj (pi) ^ V2mjv (crfc) 

V 27^77 <^ 



(159) 



ll>Nj'{p'l) ^ \/2miV <PNj"- ( PNi' 



2m N 



-jerk') 



£' = 3,4: V^' (-P2) - ^V^m" 2m, Xa^' 



■C+ -(°- fe ) v 

2m^ ' 



c+ 



where (£>jvj and %a£ are two-component spinors. For the 
waves with J = 0, 1 we have 

L = 0, J = : (^b'OO 000 ^')^!-^)) 

x (^(-P 2 )0° 00 (A:)^(p 1 )) A%fJ» A {s) * 
(x5w^Wj) V4™/v™A 



,(0,00,0) , , 



L = 0, J = 1 : (^(piJQi 01 ^)^-^) 

x (^A(-p 2 )Qi 01 (fc)^( P1 )) ^SaM * 



xA 



(1,00,1) 
N A—*N A 



(S). 



(160) 



Let us consider bispinors with real components. For nu- 
cleons, we write 

= (nlm) ' ^ = ( ^ (JVAw(JVi)) • (161) 



For the A, we determine the bispinor to be given by 
V T (AM _ / ^(^) \ 



(162) 



Within this definition, we can re-write (160) in terms of 
the traditional technique which uses the Clebsch-Gordan 
coefficients. For J = 0, we have 



1 

71 



(<Pi(Nj) n (Ae)-<p l (Nj)^(Ae)) = 



and for J = 1 , J 3 = 
V2 



-+ CT 3 \ / + CT 3 

XAi-^PNj ) = ( Vhj-jgCA, 



-j= (<p t W)<Pi(M) + n(Nj)^(Ae)) = 

E C 1/V l/2-.^(%-aM. 



(164) 



(ii) The D-wave component in the operator 7^. 

Equations (159) and (160) allow one to see easily the 
existence of the D-wave admixture in the operator 7^. 

By using the operator Q^ 01 (fc) = 7^ in (160), one has the 
following ncxt-to-leading term in the (J = l)-wave: 



(crfc') (crfc') 
cr 

2 771 j\j 2mA 



Xav 



VW^^Sa(^)- (165) 



Xai 



The spin operators in (165) can be presented as 



(crfc) (crfc) fc(crfc) 

-cr- ~ — — + ctO 

2mA 2m^r 2mAmjy 



k 1 



m A mM 



, (166) 



where the first term in the r.-h. side refers to the D-wave, 
while the second one gives the correction to the S'-wave 
term. In the operator r a (kj_), the D-wave admixture is 
canceled due to the second term: — [4sfc^ Q (fc^7)]/ [(itin + 
m A )(y/s + m N + m A )(s - (m N - m A ) 2 )} 



Appendix 3. Useful relations for „ and 

lS2...V n 

In this appendix, we list a few useful expressions. 

In- 1 



ya _ Y'(n— 1) f)Ccv 2 



a^k 1 - 



rt 1 - P' 



9a 9/3 



-fcT 



-fn-2 2P n _ 1 ) + 



-,P" 



a 

(n+l) 



(167) 



P" 



(168) 



x I ^l k l\/l± 



n+l 



n 1 - P' 



+ 



+ 



9" *73 f p" _ 9 p/ \ 



-,p" 



z a , 1 ...,M(- l ) n o^:x fivl ... Vn {k) 



OL n -\ 

n(n + 1) 



pit 

(169) 
(-1)" 



x(-fci)h/ fc iV9i 



n+l 



A 



P" 



"-'n+l mr- j—rt-Ln ' r^r r . 



(170) 



We now consider some further expressions used in the 
one- loop diagram calculations. In our case, the operators 
are constructed of X^ 1 }^ and Z@ im _ , where a and [3 
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indices to be convoluted with tensors. Let us start with 
the loop diagram with a Z-operator: 



For different tensors T a p, one has the following A's: 



T a p — g a 0, 

T a /3 — k^kp, 



n 



n i,„ i2n-2 



A = 



a. 



2n+V 



(172) 
(173) 



Equation (172) can be easily obtained using eqs. (167) and 
(50), while Eq.(173) can be obtained using eqs. (42) and 
(50). For the X operators, one has 



/ 



dQ 

47T 



X&+% n {^)T aP X<£X = \0%;;;? : (-1)", 

(174) 



where 



Ta/3 — 9a/3, 



A = -^^|fc| 2 "+ 2 , 
n + 1 



Tap — k^kp , 



r^| fc | 2n +4 



2n + V 



To derive Eq. (174), the properties 



aVl - Vn 2n+l 



(175) 



(176) 



of the projection operator and Eq. (43) are used. The in- 
terference term between X and Z operators is given by 



df2 

47T ' 



(177) 



with 



T Q/ 3 — <? Q( 3, 
7a/3 = fc^fc^, 



A = 0, 



^| fe |2„+2_ (1?8) 

2n + 1 1 1 v ; 



Eq. (178) is calculated using Eq. (167) and the orthogo- 
nality properties (44) of the X operators. 



Appendix 4. NA one-loop diagrams 

The calculation of the one-loop diagram for different ver- 
tex operators is an important step in the construction of 
the unitary NA amplitude. Consider the loop diagram for 
S = 1 and derive all expressions in case of different parti- 
cle masses (mi is mass of A and m 2 is nucleon mass). 

Let us start with 3 Lj (J = L — 1) states. Using the 
expression 



Sp i^s (mi +fci)(^ fel 
K 9a 



)«75(m 2 -fc 2 ) 
)(s-S 2 ), (179) 



where 5 = mi — m 2 , the one- loop diagram for the operator 
(100) is given by 



/^^i...M»( fc± ) T ^<..-»( fe± ) = A ^X B (- 1 ) n -( 1 71) | 



47T 



^Jrni + fci)^ 



; (S 6 j nTTV i+ m 2 (2n+l)J |fe| 



3 (S ^ j n+ll i+ m 2 (2n+l), 

x« (-I)"- 

Here, eqs. (174) and (175) were used. 
For 3 Lj (J = L + 1) states, one has 



(180) 



dfl 
Air 



Sp 



<K ( 1 2) .l(m 2 -fc2)] 
W 2 



_Lfci -Lfei 



) 



3 (S 6 } 2n-lV + ml(2n+l)) m 

xO^;X" (-!)"■ ( 181 ) 

Direct calculations also show that transition loop dia- 
grams between 3 L,j (J = L — 1) and 3 Lj (J = L + 1) 
states are equal to 



/ 



-Lfci^-Lfci 



4 a i lfrl 2 "+ 2 
_( 8 S 2 ) 11 ''"<- '>" + -' 

3 2/Z "I - 1 772-^ 



-0£.X(-i) r 



(182) 



One can also introduce the pure spin operator in a way 
that the transition loop diagram is equal to zero. Then 
eqs. (100-102) can be rewritten in the following way: 

W^ n = Mki)r% 2 V^ n u(-k 2 ) , i = 1, 2 , (183) 
where 



r 3/2 

L a fi =5a/3 



(.s + Af<5)(Vi + Af)(Vi + <5) 
Then, it is easy to find that 



• (184) 



Sp 



_Lfei _Lfei 

3 



xi-f 5 (m 2 - k 2 ) 



g a0 (s - S 2 ) . 



(185) 



Thus, the transition loop diagram vanishes identically. 
For 3 Lj (J = L) states, one has: 



xV^i n (m 2 -k 2 ) 



4 77 4- 1 

-(« - ^ 2 ) sa „_ li ^|fc| 2 "0^;(-l)". (186) 
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To calculate loop diagrams with 5 = 2, the following ex- 
pression is used: 



_l_ic -Lfei 

, fcl _ _ 

J a o 



xiv 2 (m 2 - k 2 ) 
a.O^l + a 2 Zi ia2 Zl i02 + a 3 Xi% 2 Xj» 02 , (187) 



U 0102 



where 



a x = 2(s - 5 2 ) , a 2 = 



326 



16 



9mi 27mf 



a 3 



64 
27m? 



(s - (mi + m 2 ) 2 ) , 
(188) 



For 5 ij ( J = L + 2), the operator one-loop diagram is 
equal to 



/ 



dQ 

47T 



Sp 



7^7^ 



| fe( ^-4 ( _ iro „ 



2n - 3 

x ( Oi + - 
V 4 2n- l v 11 2n+l 

For 5 Lj (J = L — 2), the one- loop operator is given by 



" '" ^(- a2 |fc| 2 + a3 -^Tl fc | 4 ))- (189) 



/f^K 5) ^(rni+^)(C 

x^.K^-fe)" = 

a „|fc| 2 "+ 4 (-l)"0^;;;^ 



7^7^ 



(2n + 3)ai 9/ a 2 |fc| 2 a 3 |fc| 



+ 



+ 



+ l)(n + 2) 4 V n + 1 ' 2n + 
while for 5 Lj (J = L), the operator is written as 

r dQ 



J-)), (190) 



/ 



4tt 



Sp 



±k, J-fci, 

_Lfei '/3 



«n-l 



|fc| 2 "(-l)"0"; 



2n(2n+l) 1 '" 1 v ~' 1/1 

■(2n + 3)(n+ l)ai 9 

. 3^ 8 

, a 3 |fc| 4 (n+l) 2 ^ 



|fc| 2 a 2 ( 



/2n + 5 2n+l 



+ 



2(2n- 1) 



n(2n-\)) 
(191) 



The one-loop transition diagram between 5 Lj (J = L + 2) 
and 5 X,/ ( J = L — 2) states can be expressed as 



-Lfci^-Lfci 



/ d ^S P [v^ n ( mi + k)(C> - (192) 

x< 5) l(™ 2 - fc 2 )] = 9 -^a 3 \k\^H-iro^: 



and the one-loop transition diagram between 5 Lj (J 
L + 2) and 5 Lj (J — L) states as 



/ 



dQ 



xV^i n (m 2 -k 2 ) = 
:},!; ' ::{n |fc| 2 "(-l)"0^;;.f 



±fci -Lfei 
la If3 



8(2n + l)(2n- 1) 

/2n + 3 

x a 2 

\ n 



2|fc| 2 a 3 ) 



(193) 



For the one-loop transition diagram between 5 Lj (J = 
L-2) and 5 Lj (J = L) we get 



dQ 

47T 



Sp 



xV^i(m 2 -k 2 ) 



_Lfci J-fci. 
_Lfci la 10 



3a r < 



\2n+i 



(194) 



8(2n+ 1)' 
and, for 5 ij (J = i — 1), we find to 



/ 



47T 



±fei /v -Lfei 



sa„_i 



2(2n+ 1) 



xV^i(m 2 -k 2 ) 



| fc |2n+2 ( _ 1) n Mi ; ..M» 



oi (n+ l)(2n 2 +n- 2) 9 2 n+1 
n 2 (2n-l) 8 1 1 fl2 2^I 



(195) 



Finally, the operator one-loop diagram For 5 Lj (J = L+l) 
is equal to 



/ 



dQ 
4n 



Sp 



x< 8) l(m 2 -fc 2 )" = 
sa„_ 2 (n + 1) 



7^7^ 



2{2n- l)(2n-3)) 
X ( ai -^ |fc|2fl2 2nTl) 



|fe| 2 ™- 2 (-l)™0^;;;/f" 



(196) 



and the one-loop transition diagram between 5 Lj (J 
L — 1) and 5 ij (J = L + 1) can be written as 



dQ 



xV^i n (m 2 -k 2 ) 



T- Lfcl 'Y7 fel 
J_fcj _ fa f/3 

3 



(197) 



4n 2 



^|fe| 2 "(-l)"0-;.V- + ^|fc| 2 o 2 (n + l)) . 
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Appendix 5. Amplitude of the triangle diagram 



• ► P 4 




P-Pa 



Fig. 7. Triangle diagram. 



In the last two appendices, we give results on trian- 
gle diagrams in numerical form. First, we calculate the 
triangle-diagram integral which enters Eq. (126): 



A ml-kl-i{) 



(198) 



TO^-(p-p zi + £; 7r ) 2 - im A F A m 2 N - (pa - K) 2 - «0 ' 
Notations of the momenta are illustrated by Fig. 7. Here, 
p = Pl +p 2 , p 2 = W 2 , p 2 A = s. (199) 
The physical region is located in the interval: 

(m N +m n ) 2 < s < (W-m N f. (200) 

The triangle-diagram amplitude ^triangle > s ) de- 
tcrmincd by (198) is shown in the physical region (200) in 
Fig. 8 (left column). In the right column, there are posi- 
tions of the logarithmic singularities on the second sheet 
of the complex-.s plane. Physical region of the reaction is 
also shown (thick solid line): it is located on the lower edge 
of the cut related to the threshold singularity (thin solid 
line). The positions of logarithmic singularities read: 

>)_ m 2, ra 2 (W 2 -M\ - m 2 N ){Ml +ml- m%) 



(±) 



2M% 



± 



(ml - (M A - m N ) 2 )(ml - (M A + m N ) 



1/2 



x(W 2 - (M A -m N y){W 2 - {M A + m N Y) . (201) 

Here M\ = m A — im A r A . 

In the left column of Fig. 8, the real and imaginary 
parts of the amplitude (198) at different total energies 
W are shown by solid and dashed curves, respectively. 
In the right column, one sees the singularity positions, 
s f-) (black circles) and (black squares). When 
dives into the third sheet, its position is shown as an open 
square. 



Appendix 6. Amplitude of the box diagram 
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Fig. 8. Triangle diagram amplitude. In the left columns, real 
and imaginary parts of the amplitude are shown by solid and 
dashed curves, correspondingly. Initial energy, W, is shown 
on the top of each panel. In the right columns, singularity 
positions, Eq.(201), are shown on the 2nd sheet of the 

complex-s plane. When dives to the 3rd sheet, its posi- 
tion is shown by open square. 



Here, we calculate the box-diagram integral which enters 
Eq. (136), the notations of momenta are given in Fig. 9. 
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Recall that in (205) and (206) M\ is given by 
M 2 A = m 2 A - iniA^A- 
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(207) 



Fig. 9. Box diagram. 



, S3, s 4, Sl2) = 

d A k' 1 



i(2tt) 4 (m\ - {\p + k'+ p 3 ) 2 - tm A r A ) 
1 

(m A - (\p - k' + p A ) 2 - imA^A) 



(m 2 - (\p + k'Y - i0)(m 2 - {\p - k'Y - iO) 



. (202) 



Remind that s 3 = (p-p 3 ) 2 , s 4 = {p-pi) 2 , s 12 = {P1+P2) 2 , 
W 2 =p 2 . 

In Fig. 10, we show the results of our calculation of 
A^™ lcss (W 2 , S3, S4, S12) as a function of pion-pion energy 
squared si 2 at different total energies W, under the fol- 
lowing constraint on S3 and s 4 : 



S3 = S 4 = V*12 W + m N ■ 



(203) 



This constraint corresponds to the following kinemat- 
ics in the cm. system: 



P 
Pi 

P2 
P3 

Pa 



^;0;0;0), 



pL;0;0;pi z ), 



\Jml + p 2 lz ;0;0;-pi z ) 1 
\Jm 2 N +pj z ;0;0;p 3z ), 

y tr n 2 N + PL;0;0;-p 3z ), 

^/m 2 7r +p 2 z + ^ m 2 N +p 2 z = W/2 (204) 

The positions of the box-diagram singularities are given 
by the formula: 



1 



s\T = 2m 2 . + t^(s 3 - m 2 N )(s 4 - m%) 



(205) 



+ 



>12 ' 214 72 * 

(2W 2 M 2 A -W 2 {s 3 - m 2 N )){2W 2 M 2 A -W 2 {s A - m%)) 



2W 2 {{W 2 - 2M 2 A ) 
r2 ji/f2 



2 - 4M 4 ) 



(s 3 -m 2 N ) 2 2 (2M/ 2 Mj-iy 2 ( S3 -m 2 f )) : 
2VK 2 71 2^ 2 ((iy 2 -2M 2 ) 2 -4Af4 

(s 4 -m 2 r ) 2 2 (2^ 2 M 2 -^ 2 ( S4 -n4)) 2 
2VK 2 "** 2W^ 2 ((iy 2 -2M 2 ) 2 -4M4) 

At s 3 = S4, Eq. (205) reads 

^ 2 (2M 2 -sg + m 2 ,) 2 



box 
*12 



: 4m, 



(W 2 - 2M 2 A ) 2 - 4M% ■ 



(206) 
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Fig. 10. Box diagram amplitude as a function of S12 under the 
constraint (203) (corresponding magnitudes of S3 and S4 are 
shown in the right column. In the left columns, real and imag- 
inary parts of the amplitude are shown by solid and dashed 
curves correspondingly. Initial energy, W , is shown on the top 
of each panel. On the right columns singularity positions, s\%*, 
Eq.(205), are shown on the 2nd sheet of the complex-si2 plane. 
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